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Abstract 

The Li(p, 7o) 8 Be reaction is studied in a microscopic cluster model. All rel- 
evant subsystem properties are well reproduced. The calculated astrophysical 
S factor is in good agreement with the experimental data, although some Ml 
strength is missing in the 1 + ; 1 resonance region. There is no contribution 
from the 1 + ; state. We estimate an Ml strength of 3.5% at 80 keV. By as- 
suming the experimental resonance parameters and channel spin ratio for the 
first 1 + state that Ml strength is 6.3%, in good agreement with an R matrix 
analysis. The large (20%) Ml strength determined from a transition matrix 
element analysis may be caused by assuming a vanishing OP3/2 spectroscopic 
factor, which is not supported by shell model calculations. 
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I. INTRODUCTION 



Astrophysical processes most often involve charged particle nuclear reactions at very low 
energies. The cross sections of such reactions decrease exponentially with decreasing energy, 
making their laboratory measurements extremely difficult. One of the most important reac- 
tions in nuclear astrophysics is the radiative proton capture on 7 Be. This reaction produces 
8 B in the Sun, whose (3 + decay is the main source of high-energy solar neutrinos. That 7 Be 
is radioactive makes the measurement of the low-energy 7 Be(p, 7) 8 B cross section even more 
difficult. There is a considerable spread among both the experimental (see, e.g., 0) and 
theoretical (see, e.g., [0]) values of the low-energy astrophysical Syj factor. However, a com- 
mon feature in all of these analyses is the assumption that the cross section at low energies 
(E < 200 keV) is dominated by the El transition. Hence the extrapolated astrophysical S 
factor is flat. 

The 7 Li(p, 7) 8 Be reaction was studied for energies from to 80 keV by Chasteler et al. ]3[] 
with the hope that the somewhat easier experimental situation ( 7 Li is stable) might allow one 
to gain some insight into the similar 7 Be(p, 7) 8 B process. A considerable forward-backward 
asymmetry was found in the angular distribution of the emitted photons, indicating a non- 
negligible Ml (p-wave) strength at such low energies. A transition-matrix element (TME) 
analysis of the data was performed, resulting in two possible values of the Ml strength 
of 20% and 80%. It was suggested that the situation may be similar in the case of the 
7 Be(p, 7) 8 B reaction. If that were so, then the lowest energy (^ 100 keV) experimental 
cross section values would contain a 20% (80%) Ml contribution. Thus, as the Ml cross 
section becomes negligible at lower energies because of the centrifugal barrier, the 5* factor 
would not be flat but decreasing with decreasing energy, and the zero energy Su factor 
would be 20% (80%) smaller than currently believed. 

Rolfs and Kavanagh j|] suggested that the observed low-energy Ml strength of the 
7 Li(p, 7) 8 Be reaction comes from the tail of the 1 + state at Eiab = 441 keV. More recently, 
the astrophysical 5* factor of the 7 Li(p, 7) 8 Be reaction was accurately measured between 100 
and 1500 keV by Zahnow et al. [||]. The data were interpreted again by assuming that the 
small Ml contribution comes from the tail of the 1 + resonance. Weller and Chasteler |5| 
pointed out that some of the formulae used in the previous analyses |||| were incorrect, 
and that the resonance tail could only result in a 2% p-wave contribution, much smaller 
compared to the earlier estimates. Barker ]7| achieved acceptable fits to the measured 
angular distributions and analyzing powers by attributing the p-wave Ml contribution to 
the tails of the low-energy resonances. He obtained an Ml strength of 5 — 6%, again much 
less than that reported by Chasteler et al. However, to achieve these fits, some of the 
spectroscopic amplitudes needed the opposite signs to those obtained from the shell model. 

Barker also noted J7| that the analogue of the 7 Be(jo, 7) 8 B reaction in 8 Be would be 
the one which leads to the 2 + ; 1 isobaric analogue state of the ground state of 8 B. Such a 
measurement was performed by Godwin et al. || and no indication of any low-energy Ml 
strength was found. Therefore, it may be assumed that the standard extrapolation of the 
7 Be(p, 7) 8 B S factor is correct and reliable. However, the origin of the discrepancies among 
the various analyses of the 7 Li(p, 7) 8 Be reaction is still not understood. 

Herein, we study the role of the Ml transition in the 7 Li(p, 7) 8 Be reaction by using a 
microscopic eight-body model. 
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II. MODEL 



Our model is a microscopic multi-cluster ( 4 He+ 3 H+p; 4 He+ 3 He+n; a+a) Resonating 
Group Method (RGM) approach to the eight-nucleon system. The trial function of the 
eight-body system is 



*= e ^{[[$ p (^)] s XN 2] L(Pi,p 2 y 

S,li,h,L 

+ E A{[[* n (* a &)] a xfrw(jh,t>*) 

S,h,h,L 

+ A{[^xl(p)} jm }, 



JM 



JM 



(1) 



where A is the intercluster antisymmetrizer, the $ cluster internal states are translationally 
invariant harmonic oscillator shell model states (a = 4 He, t = 3 H, and h = 3 He), the p 
vectors are the different intercluster Jacobi coordinates, l\ and l 2 are the angular momenta 
of the two relative motions, L and S are the total orbital angular momentum and spin, 
respectively, and [...] denotes angular momentum coupling. The sum over S,h,l 2 , and L 
includes all physically relevant angular momentum configurations. In order to project out 
pure 7 Li and 7 Be states, the angular momenta are recoupled in the 7+1 channels following 
the scheme 



[(h, l 2 )L, (si, s 2 )S]J -> [(h, si)I 7 , s 2 ]I, h 



J. 



(2) 



where is the total spin of the seven nucleon system, and I is the channel spin. 

The 7 Li(p, 7o) 8 Be radiative capture cross section is calculated perturbatively 0, and is 
given by 
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where Ij and s 2 are the spins of the colliding clusters, A denotes the rank of the electro- 
magnetic operator Ai^ (fl = E or M), u represents the entrance channel, and Jf and J, is 
the total spin of the final and initial state, respectively. The initial wave function ^ t * T is a 
partial wave of a unit-flux scattering wave function. We follow the convention of Edmonds 
for the reduced matrix elements. 

The final 8 Be ground state is unbound in the a + a channel which, in principle, could 
cause convergence problems. However, this is not the case here since the 7 Li+p and 7 Be+n 
channel wave functions are orthogonal to the a + a one in the asymptotic spatial region. 
This feature allows us to treat the 8 Be ground state as a pseudo-bound state in the a + a 
configuration. The orthogonality also indicates that the 7 Li+p — > a + a contribution to the 
cross section is small, and that the inclusion of the 7+1 configurations in the 8 Be ground 
state is important. 

Putting ([I]) into the eight-nucleon Schrodinger equation which contains the nucleon- 
nucleon (NN) strong and Coulomb interactions, we arrive at an equation for the intercluster 
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relative motion functions x- F° r the + bound state these functions are expanded in terms 
of products of tempered Gaussian functions exp(— 7«p 2 ) [fl ij with different ranges ji for each 
type of relative coordinate. The expansion coefficients are determined from a variational 
principle. The scattering states are calculated from a Kohn-Hulthen variational method 
|TT| for the S'-matrix, which uses square integrable basis functions matched with the correct 
scattering asymptotics. Then, using the resulting eight-nucleon wave functions, the cross 
section (ED is evaluated. 



III. RESULTS 

The model used herein is essentially the same as that used in Refs. [|T^,|T^] , which was 
used to describe the 7 Li, 7 Be, 8 Li, and 8 B nuclei, and the 7 Be(p, 7) 8 B reaction. In those 
calculations all possible arrangements of the three clusters were included, e.g., (ah)p, (ap)h, 
and (hp)a for 8 B. However, we consider only the 7 Li+p and 7 Be+n configurations due to 
computational constraints. The results of a full model and a model which contains only 
7 + 1 configurations are quite similar |Tj| , provided the experimental subsystem properties 



(separation energies, channel thresholds, sizes, etc.) are correctly reproduced in each. 

We use the same parameters for the description of 7 Li and 7 Be, and the same NN 
interaction (Minnesota force) as in Ref . |12| . Using these parameters, the 7 Li and 7 Be prop- 
erties are well reproduced fl2|l . The calculated relative binding energies 7 Li(g.s.)- 7 Be(g.s.) 
and 7 Li(g.s.)- 8 Be(g.s.) are also close to the experimental values. Note that in eq. (Q) the 
experimental value for E 1 is used. 

We calculate the El, E2, and Ml cross sections for the capture to the ground state of 
8 Be. An unambiguous treatment of the transition to the broad 2 + excited state of 8 Be would 
require a more sophisticated model. 

As the final state is J n = + , the El, E2, and Ml transitions occur from 1~, 2 + , and 1 + 
initial scattering states, respectively. The 1~ and 2 + partial waves are non-resonant above 
the 7 Li+p threshold, while there are two 1 + resonances at E c m = 386 keV (r = 11 keV, 
T = 1) and E c m _ = 897 keV (r = 138 keV, T = 0). In order to explore the role of the 
Ml capture at low energies, the properties of the 386 keV state must be reproduced. We 
decreased the exchange mixture parameter of the NN interaction by 0.4% and achieved 386 
keV position and 13 keV width. The energy and width of the second 1 + (T = 0) state 
are 1050 keV and 342 keV, respectively, in our model. These resonance parameters were 
determined from an analytic continuation of the scattering matrix to complex energies |T4 



As the T = state is higher in energy than experimentally observed, and much broader, it 
does not have any observable effect on the low-energy cross section. This is not a serious 
problem, as only the first state is expected to contribute significantly to the cross section at 
~ 80 keV @||. Note that the 0.4% change in the NN interaction has a negligible effect on 
the subsystem properties. 

We generate the + bound state and 1~, 2 + , and 1 + scattering states using the same force 
in all cases. In order to account for specific distortions of 7 Li and 7 Be, we describe these 
nuclei by using six basis states. The resulting large number (50-60) of coupled channels 
required special numerical attention to keep high precision in the scattering calculations. 
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The result of the calculation of the 7 Li(p, 7) 8 Be cross section is presented in Fig. [I] in 
terms of the astrophysical S factor 

S(E) = a(E)E exp \2irr){E)\ , (4) 

where rj is the Sommerfeld parameter. As our E2 cross section is 2-3 orders of magnitude 
smaller than the El, it is not displayed. There is no contribution to the cross section coming 
from the second 1 + state, as expected. The nonresonant part of the cross section and the 
first 1 + state is well reproduced, although some Ml strength is missing in the resonance 
region. 

Note that, as discussed above, inclusion of the 7 Li+p and 7 Be+n configurations in the 
ground state of 8 Be is necessary. The contribution to the cross section from the a + a 
final state configuration is very small. This + state is thought to be one of the pure two- 
cluster (a + a) systems JTH] . However, by adding the high- lying 7+1 channels to the a + a 



configuration, one gains 1-2 MeV of energy in the + state, indicating that this state may 
not be a pure a + a configuration. 

In our model the 386 keV (T = 1) resonance is dominated by the 1 = 1 channel spin. 
This is in contradiction to experiment fllBl , which gives 3-5 for the ratio of the / = 2 to 
1 = 1 channel spin contributions to the cross section. The same 1=1 dominance was 
also observed in the calculations of the 1 + isobaric analogue state in 8 Li and 8 B. In order 
to understand the origin of this 1 = 1 dominance, we calculated the percentages of the 
orthogonal {LS) components in 8 Li (which is the only two-body bound state in the triplet). 
The (LS) = (1,0), (0, 1), (1, 1), and (2, 1) components have 83.4%, 0.03%, 13.3%, and 3.2% 
weights, respectively. However, the contribution of the largest component, (1,0), to the 
1 = 2 channel spin state is zero, and hence the 1 = 1 configuration dominates. Very similar 



(L, S) components were found in a different three-body model |17| using completely different 
cluster-cluster interactions. That model would also predict a dominant 1 = 1 component in 
the A = 8 1 + ; 1 triplet. In the second 1 + (T = 0) state, our model predicts a dominant 1 = 2 
component. It remains a question as to whether the mixing of these states could influence 
the experimental determination of the channel spin ratio. 

As a comparison to the results obtained from our RGM model, the I = 1, 2 spectroscopic 
amplitudes were calculated using wave functions obtained from a large space shell model 
calculation. The interaction used was that of Zheng et al. ||18|| , based on a multi- valued 
G-matrix calculation and defined within the complete (0 + 2 + 4)huj shell model space. The 
calculations were performed using the shell model code OXBASH ||19|| . The spectroscopic 
amplitudes Sj=i,2 are defined in terms of the shell model spectroscopic amplitudes by 

S I = J2 a i,j$ji (5) 

j 

where dj are the shell model spectroscopic amplitudes (j = 3/2 or 1/2), 

0j= ,] (j f ( s Be)\\a]\\l 7 fU)), (6) 
' 2J f + 1 



and aij are transformation coefficients for the recoupling of the angular momenta in terms 
of the channel spin, /, 
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J f . (7) 



Table | lists the spectroscopic amplitudes, Sj, as obtained from our shell model calculation, 
along with those obtained from the calculation using the Cohen and Kurath (CK) (6 — 
16)2BME Op-shell interaction [KJ. Note that the contribution from higher shells to the 
amplitudes calculated in the (0 + 2 + A)tiu space is negligible. The CK amplitudes favor 
the 1 = 2 channel spin for the 1 + ; 1 state, while the Zheng ones support the results of our 
analysis. This discrepancy is the result of a slight redistribution of particle strength between 
the 0p 3 /2 and Opi/2 orbits, as is illustrated in Table |TJ. The spectroscopic amplitude, 9j, for 
the capture to the Op^/2 orbit is larger in the Zheng model, while there is a corresponding 
decrease in the magnitude for the capture to the 0pi/ 2 orbit. It is also interesting to note 
the signs of the channel-spin amplitudes for the 1 + ; state. In Barker's two- level i?-matrix 
calculation 0, the signs of the reduced widths used were obtained from the shell model. 
However, Barker changed the sign of the 1 = 1 amplitude relative to that of the 1 = 2 one, 
in order to fit the data. Our shell model results does not support this sign change. 

We note, that if our channel spin ratio is in fact wrong, then probably the singlet-odd 
component of the NN force is responsible for this failure. In our model we calculate the 
peak to 80 keV S factor ratio for the 1=1 and 1 = 2 channel spin contributions. Using 
these numbers we can estimate the Ml strength at 80 keV for any given channel spin ratio. 

The necessity of the microscopic description of the 7 Li(p, 7) 8 Be (in contrast to that for 



Be(p, 7) 8 B [21], for example) is demonstrated by noting that in a potential model the 1 = 2 
contribution would come almost exclusively from the internal magnetic moments of 7 Li and 
7 Be. The orbital part of the Ml operator describing the relative motion (which strongly 
dominates the spin part) is a rank-zero tensor in spin space. Thus, for any non-vanishing 
cross section from the 7 Li-p relative motion, the initial and final channel spins must be equal. 
Hence, as the channel spin in the + is / = 1, the 1 = 2 orbital cross section is zero. In our 
microscopic model, the 1 = 2 cross section comes mainly from nucleon exchange effects. 

The S factor is displayed on a linear scale in Fig. [| The nonresonant El component is 
in good agreement with the data. We emphasize that in our model all parameters are fixed. 
The requirement that all relevant properties of the subsystems be reproduced does not leave 
any room for adjustment of parameters. 

The contribution of the Ml component to the total S factor at 80 keV is 3.5%. If one 
adopts a value of 3.2 for the 1 = 2/1 = 1 channel spin ratio, together with the experimental 
width and peak value of the 386 keV resonance, this Ml contribution is 6.3%, consistent 
with the estimates by Barker f7|. The value of 20% determined from the TME analysis 
may be due to the assumption of a negligible 0p 3 / 2 spectroscopic factor. That assumption 
is not supported by the results of the shell model calculations presented in Table [H|. 



IV. CONCLUSIONS 

In summary, we have studied the low-energy 7 Li(p, 7 ) 8 Be reaction in a microscopic 
cluster model. All relevant subsystem properties are well reproduced by our model. The 
parameters of the 1 + ; T = 1 resonance in 8 Be are in good agreement with the experimental 
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values, while the 1 + ;T = state is much higher and much broader in energy compared to 
experiment. 

We have calculated the El, E2, and Ml cross sections perturbatively. The nonreso- 
nant El component agrees well with the experimental data, while the E2 contribution is 
negligible. The calculated full astrophysical S factor is in good agreement with the data of 
Zahnow et al. ||, except that the l + ;0 state does not contribute in our model, and some 
Ml strength is missing in the 1 + ; 1 resonant region. This is partly caused by the fact that 
the width of the 1 + ; 1 state is slightly larger compared to experiment. 

Our model predicts that in the A = 8, 1 + ; 1 triplet the dominant configuration has 
1=1 channel spin. Consequently, the 7 Li(p, 7 ) 8 Be Ml cross section is also dominated 
by the 1 = 1 component at E cm . = 386 keV, in contrast with the 1 = 2 dominance 
observed experimentally. Results of large-space shell model calculations, which give a larger 
spectroscopic factor for the 1 = 1 component in the 1 + ; T = 1 state, seem to support our 
RGM result, although this depends sensitively on the distribution of strength between the 
0p 3 /2 and Opi/2 orbits. We note, however, that the total Ml cross section can be estimated for 
any channel spin ratio in our model, using the calculated 1 = 1 and 1 = 2 Ml components. 

The calculated Ml strength at 80 keV is 3.5% in our model. If we assume the experi- 
mental peak height and width of the 386 keV resonance in the S factor, and the 3.2 1=2/1=1 
channel spin ratio, the model predicts 6.3%, in good agreement the result of an R matrix 
analysis @. 

Our shell model calculations show that the OP3/2 7 Li+p spectroscopic factor is larger 
than the Opi/2 one in the 1 + ; 1 state. Thus the assumption of a vanishing spectroscopic 
factor in the TME analysis is incorrect. This may have led to the large Ml strength found 
therein. 
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FIGURES 

FIG. 1. 7 Li(p, 7o) 8 Be astrophysical S factor as a function of the 7 Li+p center-of-mass energy. 
The result of the present calculation is given by the solid line, while the El and Ml contributions 
are given by the dashed and dot-dashed lines, respectively. The data are those of Zahnow et al. 

FIG. 2. As for Fig. |, but displayed on a linear scale. 
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TABLES 



TABLE I. Spectroscopic amplitudes, Si, labeled by the channel spin / for the 7 Li(p, 7) 8 Be 
reaction to the 1 + states in 8 Be, as obtained from the shell model calculations using the interactions 
of Zheng et al. ]l8| and of Cohen and Kurath pCfl . 



State 


E c . m , (keV) 


I 


Si 

Zheng 


CK 


i + ;i 


386 


1 


0.3499 


0.2644 






2 


0.2774 


0.3910 


l+;0 


897 


1 


-0.2612 


-0.2773 






2 


0.5132 


0.5623 



TABLE II. Spectroscopic amplitudes, 9j, for the 7 Li(p, 7) 8 Be reaction to the 1 + states in 8 Be, 
as obtained from our shell model calculations. 



State 


E c . m . (keV) 


3 


Zheng 




CK 


i + ;i 


386 


0^3/2 


0.4327 




0.4009 






Opi/2 


-0.1104 




-0.2490 


l+;0 


897 


0P3/2 


-0.0290 




-0.0334 






Opi/2 


-0.5752 




-0.6266 
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